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Recent developments in lattice QCD calculation of weak matrix elements involving light quarks are described. 
We focus on four topics: Bk with the Kogut-Susskind and Wilson quark actions, A/ — 1/2 rule, proton decay 
matrix elements and the application of domain wall QCD to calculation of weak matrix elements. 
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1. Introduction 

The kaon bag parameter Bk, which is re- 
quired for constraining the CP violation param- 
eter in the Cabibbo-Kobayashi-Maskawa matrix 
from experiment, has been successfully calculated 
with the Kogut-Susskind(KS) quark action, tak- 
ing advantage of its Ua(1) chiral symmctry[|l|. 
A recent systematic calculation of Bk with high 
statistics Q, however, reveals that the system- 
atic uncertainties in the perturbative renormal- 
ization factors to connect lattice and continuum 
operators have sizable magnitude, which indi- 
cates the necessity of non-perturbative renormal- 
ization. For the calculation of Bk with the Wil- 
son quark action we have to deal with the opera- 
tor mixing problem caused by the explicit chiral 
symmetry breaking in the action, which turned 
out to be not effectively treated by perturba- 
tion theory I3I. Several non-perturbative methods 
have been proposed to control the operator mix- 
ing. 

The CP conserving K — > tttt decay amplitudes 
give an intriguing testing ground for lattice QCD 
calculations. It would be exciting to be able to 
show that QCD explains the A/ = 1/2 rule: 
A0/A2 « 22 where ^0,2 = A{K -^ t:t:[I = 0, 2]). 
Although studies oi K ^ tttt decay amplitudes 
with lattice QCD were initiated a decade ago, lit- 
tle progress has been achieved since then. Re- 
cent work on the one-loop calculation of K — > tttt 
decay amplitudes using chiral perturbation the- 
ory (ChPT) |0,H|, however, provides some impor- 



tant information about the lattice results. 

In lattice QCD studies for the proton decay am- 



plitude 



pioneering work estimated it 
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from the matrix element {0\O-^\p) with the aid of 
the chiral lagrangian|^,0, which were followed by 
the direct measurement of (7r°|C'^|p)Q. The re- 
sults showed an unexpectedly large discrepancy 
between these two methods: the direct method 
yielded a two or three times smaller value than 
the indirect one. This peculiar feature was con- 
firmed by JLQCD last yearH. At this conference, 
however, JLQCD pOJ pointed out that the calcu- 
lational method in refs.pJ9| was wrong, and pre- 
sented the correct estimate of the proton decay 
amplitudes with the direct method. 

The domain wall quark formulation in lattice 
QCD||ri|,^ potentially has superior features over 
the Wilson and KS quark actions: realization of 
the chiral limit on arbitrary number of flavors 
without fine tuning. Recent simulation results 
for Bk show a good chiral behavior [|l3[, which 
encourages us to apply the domain wall quark to 
the calculation of weak matrix elements. Some 
progress has been made this year numerically and 
theoretically. 

This article is organized as follows. In Sec. 2 
we present the Bk results with the KS and Wil- 
son quark actions. From the examination of the 
systematic errors in the KS result we explain 
the reason why the non-perturbative renormal- 
ization is necessary. The present status of the 
non-perturbative renormalization method is also 
discussed. Results for the K ^ tttt decay am- 
plitudes are presented in Sec. 3, where we focus 
on the one-loop effects of ChPT. In Sec. 4 we de- 
scribe an essential advance on the calculation of 



the proton decay matrix elements. This year's 
progress on the calculation of weak matrix ele- 
ments with the domain wall QCD is explained 
in Sec. 5. Our conclusions are summarized in 
Sec. 6. Results for structure functions have been 
reviewed by Petronzio [|l4| . 

2. Bk 

2.1. Bk with the KS quark action 

The K meson B parameter is defined by 

„ _ {k°\s^,{l^j5)d-sj^{l-j5)d\K°) ,. 
"" f(XO|s7^75d|0)(0|s7^75d|i^"> ' ^ ^ 

for which we quote results in the naive dimen- 
sional regularization(NDR) scheme at the renor- 
malization scale 2GeV. In the ratio of eq.(P the 
fluctuations are largely canceled between the nu- 
merator and the denominator, which enable us 
to easily achieve high statistical accuracy. Tak- 
ing advantage of this virtue, we have had control 
over the systematic errors step by step. As of this 
conference Bk is the weak matrix element with 
which we have had the most success. 

Figure |l| presents a recent result for a sys- 
tematic quenched Bk calculation performed by 
JLQCD employing seven values of lattice spac- 
ing with rud = 'TisP]. The non-local weak four- 
fcrmi operator is constructed either in gauge in- 
variant way with insertion of the link variables 
(Cinv) or in non-invariant way without link vari- 
ables in Landau gauge (Onon-inv)- Each operator 
is perturbatively matched to the operator in the 
continuum NDR scheme at the scale q* = l/a em- 
ploying tadpole improvement. The value of Bk 
at /x = 2GeV is obtained via a two-loop renor- 
malization group running from /i — 1/aGeV. 

The statistical error is 0.1% at the largest lat- 
tice spacing and gradually increasing to 1.2% at 
the smallest one. Let us examine conceivable sys- 
tematic errors in order, (i) Quark mass depen- 
dence: Bk at the physical kaon point is accessi- 
ble with interpolation, not extrapolation, which 
means the interpolated value is strongly robust 
against the choice of fitting function, (ii) Finite 
size effects: Finite size studies at /3 = 6.0 and 6.4 
show that the magnitude of the size dependence 
decreases to less than 0.5% for L > 2.2fm. The 
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Figure 1. Quenched BK(NDR,2GeV) with the 
KS quark action as a function oinipa. See text for 
the solid and dashed lines. Domain wall fcrmion 
results will be explained in Sec. ^. 

spatial size is kept larger than L ~ 2.3fm at all /?. 
(iii) Scaling violation: The five data points be- 
low mpG ~ 0.6 for each operator are consistent 
with 0{a^) scaling behavior expected theoreti- 
cally |1^. A quadratic fit of the five points(dashed 
lines in Fig |l|) gives a value at the continuum 
limit Bk = 0.616(5) for Onon-mv and 0.580(5) 
for Oinv This 7a discrepancy raises the ques- 
tion of what causes this difference, (iv) O(a^) 
uncertainties: As seen in the matching procedure 
of Oinv,non-inv to the continuum operator, their 
difference ABk should be of not only 0{a^) but 
also 0{a^). Figure illustrate that the differ- 
ence is well described by the function ABk = 
bi{mpa)'^ + b2{a^a{q*))'^ below TOpO ~ 0.6, which 
suggests that a a|^ term should be incorporated 
in the fit of Bk- An attempt to simultaneously 
fit Bk for both operators is shown in Fig. ll|(solid 
lines), which yields Bk — 0.628(42) in the contin- 
uum limit(cross). This 7% error, which is roughly 
10 times larger than those with naive quadratic 
fit, reflects the magnitude of the 0{a^) uncertain- 
ties in the operator matching procedure. To elim- 
inate the 0{a^) uncertainties, it would be neces- 
sary to adopt the non-perturbative renormaliza- 
tion method, (v) Quenching effects: In 1996 the 
OSU group studied the dynamical quark effects 
on Bk employing three ensembles of gauge con- 
figurations at a-i « 2GeV: TV/ = at /3 = 6.05, 
TV/ = 2 at /3 = 5.7 and TV/ = 4 at /3 = 5.4[|6). 
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Figure 2. Difference of 5/^ (NDR, 2GeV) be- 
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They found that Bk with TV/ = 3 is 5 ±2% larger 
than the quenched value. However, it is hard to 
assume this number as the dynamical quark ef- 
fects in the continuum limit: we cannot exclude 
the possibility that the scaling violation of Bk in 
full QCD might be different from the quenched 
case. This fact compels us to perform the sys- 
tematic calculation of Bk in full QCD varying 
the lattice spacing and the dynamical quark mass, 
(vi) Degenerate quark mass ma — rud'. For Bk 
in the quenched QCD there is a divergent chiral 
logarithm as md ^ in the chiral expansion of 
ChPTJl^,^. We can investigate the TOs ^ nid 
case only in the full QCD. 

The next step toward the precise determination 
of Bk is the control of the 0(0^) uncertainties. 
We also note that the systematic full QCD cal- 
culation of Bk with nis ^ rud is now feasible in 
view of the performance of the current full QCD 
simulations [Q . 

2.2. Non-perturbative renormalization 

A naive non-perturbative renormalization 
method is to copy what is done in pertur- 
bative renormalization in a non-perturbative 
way[|l9[. The renormalization condition in the 
momentum-space subtraction scheme (MOM) 
for a certain operator is provided by impos- 
ing that the off-shell quark matrix element in 
a fixed gauge should coincide with their tree 
level value. For example, if we consider the 
bilinear operator Or — 'ipTip, the condition 
is given by Zr{go, fJ.a){q{p)\Or\q{p))\p2^f,2 = 



{q{p)\Or\q{p))tTee, where F denotes Dirac matri- 
ces and p is the off-shell quark momentum. In 
non-perturbative renormalization on the lattice, 
the matrix element {q{p)\Or\q{p)) is constructed 
with the vertex and wave functions obtained from 
the Fourier transformed quark Green functions in 
Landau gauge. This method, which we refer to as 
NPMOM method hereafter, is expected to work 
if we can find the appropriate region for the ex- 
ternal quark momentum Aqcd <^ p <^ a~^ to 
avoid the non-perturbative contributions and the 
cut-off effects simultaneously. 

The NPMOM method has been actively be- 
ing pursued for various operators in different ac- 
tions! 20 -23. The studies, however, reveal that it 
is not a trivial task to find the appropriate region 
of p: it depends on each operator strongly; more- 
over, it is hard to find it for some operators. At 
the present stage it is not clear to what extent the 
non-perturbative contaminations and the cut-off 
effects can be controlled quantitatively. Another 
concern is the issue of Gribov copies in the Lan- 
dau gauge fixing. We have no definite way to 
estimate the ambiguities induced by the choice of 
Gribov copies. 

The technical difficulties in the NPMOM 
method can be entirely overcome by another non- 
perturbative method that uses the Schrodinger 
functional(SF) as the renormalization scheme [p5|. 
The SF scheme is essentially a finite-volume 
renormalization technique. QCD is consid- 
ered in a finite space-time volume of physi- 
cal size L in all directions, where all renor- 
malized quantities are defined at scale /i — 
1/L. A change of the lattice size at fixed bare 
parameters is a change of the renormalization 
scale. For the bilinear operator -ipTtp a pos- 
sible choice of the renormalization condition is 
Zr{go,L/a) = c^/f^/fr{t)\t=L/2, where fr{t) = 
J (fy(fz{i^f{x, t)Ti,},{x, t)Cf, {y, 0)rC/(z, 0)) and 
/cxt = / d^y'd^S'd^yd^ziCfiy', L)rC}, {z', L) 
C/' (y, 0)rC/(z, 0)) with /, /' the flavor. C, ■ • ■ , C' 
denote the boundary quark fields at i = and L 
with zero spatial momentum projection. The con- 
stant c in the renormalization condition is prop- 
erly chosen such that Zr = 1 at tree-level. 

There are three main advantages in the SF 



scheme, (i) Gauge invariance: It is not necessary 
to fix the gauge, (ii) Mass independent renormal- 
ization: This scheme ahows us to perform numer- 
ical simulations at zero quark mass. We can spec- 
ify the boundary conditions such that the lowest 
eigenvalue of the free Dirac operator is lifted by a 
gap of order 1/L at vanishing quark mass|2g] . (iii) 
Non-perturbative renormalization group: The 
scale evolution of a factor of two for the renormal- 
ization constant Z{fj,), which is expressed by the 
step scaling functions a{g'^{L)) = Z{2L)/Z{L), 
can be measured by increasing the lattice size L 
to 2L. In the next step the lattice spacing a is 
enlarged by a factor of two keeping the renormal- 
ized coupling g^ fixed such that the number of 
lattice sites L/a is reduced by a factor of two. Re- 
peating this procedure we can connect the pertur- 
bative regime to the hadronic scale: Z(Lmax) = 
Z(2-'=L^ax)a(ff'(2-'=Lmax)) ' ' ' a{f {2-^ L^,^)) 
for fc > 1. We emphasize that the step scaling 
function is defined in the SF scheme and should 
be independent of the lattice regularization. 

The ALPHA Collaboration demonstrates the 
effectiveness of this method calculating the non- 
perturbative scale evolution of the renormalized 
quark mass with Nj = in the SF scheme, where 
the scale evolution from the hadronic scale to the 
perturbative regime around lOOGeV is described 
within 1% errors |2^]. Besides the quark mass the 
SF renormalization method has been also applied 
to the twist-two operator associated with the non- 
singlet parton density |28,nj] and the static- light 
axial vector current relevant to the B meson de- 
cay constant p9|. 

2.3. Bk w^ith the Wilson quark action 

We have two purposes to calculate Bk with the 
Wilson quark action. One is to demonstrate that 
the Wilson result is consistent with the KS result, 
which would give confidence to the lattice QCD 
calculation of Bk- The other is an application to 
the heavy-light meson system for which the inter- 
pretation of flavor quantum numbers with the KS 
action is problematic. 

While steady progress has been achieved so far 
with the KS quark action, studies with the Wilson 
quark action are rather stagnant. This is due to 
a non-trivial mixing between the weak four-fermi 
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Figure 3. Quenched Sa'(NDR, 2GeV) with the 
Wilson (filled symbols) and SW (open symbols) 
quark actions as a function of m^a. For the KS 
result in the continuum limit (cross) see Fig. nl 

operator Oll = s7m(1 " 75)^ • S7^(l - 75)^ and 
other four-fcrmi operators with different chirali- 
tics. One of the mixing operators is 575^ • 575^. 
Numerical studies show that the magnitude of the 
matrix element {R'^lsj^d ■ 575^1^^'') is a factor of 
10 larger than that of {K°\Oll\K'^) at the point 
of the physical kaon mass (see, e.g., refs. [ pO| , pl| ). 
This can be reasonably expected from the ratio 
of their vacuum saturation approximations: 

16(i^0|s7,,75d|0)(0|s7,,75d|i^°) 16 ~ ^^ 

with the aid of the PCAC relation. In this 
case the mixing problem is not adequately 
treated with one-loop perturbation theory, be- 
cause two-loop contributions are potentially 
large: \al,,{K°\sj,d ■ sj,d\K°)/{K°\OLL\K°)\ 
>0.1, where we assume aiatt = 0.15 — 0.2 as the 
currently accessible strong coupling constant with 
tadpole improvement. This fact compels us to 
control the operator mixing non-perturbatively. 

Three major non-perturbative methods have 
been proposed so far. The most conventional 
method is to rely on chiral perturbation the- 
ory ||3C|,|3l[. The mixing structure of the weak 
four-fermi operator is expressed by Oll = Oll + 
"^i^iZiOi, where the mixing operators are ar- 
ranged into the Fierz eigenbasis given by Oi = 
SS + TT + PP, O2 = SS - TT/3 + PP, O3 = 
VV -AA + 2{SS - PP) and Oi = VV -AA~ 
2{SS - PP), with VV = s-f^d ■ s-f^d, AA = 
sjfijbd-sjfxjbd, SS = sd-sd, PP = sj^d-sj^d and 
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Figure 4. Mixing coefficients zi, . . . , Z4 obtained 
witli (a) the NPMOM method and (b) the Ward 
identity method as a function of p^a^ around 
strange quarlc mass with the Wiison quarfc action 
at /3 = 6.3. 

TT = sCT^^d • scr^^.d/2. The authors of ref. [|o| 
assume the generai form of the chirai expansion 
up to 0{p'^) 



{K'{Pf)\dlf 



for the matrix elements of 0£" : 
|X°(p,)) ^ a + I3ml + 



-K ^ IPi ■ Pf 

+ 5i7n\ + 52m]^pi • Pf + SziPi ■ PfY , where chi- 
rai logarithms are ignored. One-loop perturba- 
tive expressions are used for the mixing coefh- 
cients zi , . . . , 2:4 in O^^ . The coefficients in the 
chirai expansion of {K'^{pf)\0^j^l^\K^{pi)) can be 
determined by changing the quark mass and the 
spatial momentum of K^ and K'^ , with which we 
can eliminate pure lattice artifacts of a, [3 and 5i 
terms. In practical applications of this method, 
however, it is difficult to fix i5i, 82 and 83 with 
good precision. The results for Bk are plotted in 
Fig. I (triangles), which seem to be reasonable es- 
timates compared with the KS result at the con- 
tinuum limit. Nonetheless, this method is not 
promising for controlling systematic errors, since 
it contains unknown uncertainties from higher or- 
der effects of ChPT which survive even in the 
continuum limit. Another defect of this method 
is that it cannot be applied to the heavy-light 
meson system. 

The second method is to determine the mix- 
ing coefficients 21 , . . . , Z4 using the NPMOM 



method 12^. We require the Fourier transformed 
vertex function should be proportional to the 
tree-level Dirac structure: G~^{p)G^^{p) 
{dLL[ssdd]{p))G-/{p)Gf{p) _cx 7m(1 - 75) 
® 7p(l — 75), where [ssdd]{p) represents 
s{p)s[p)d{p)d[p). Ggj;{p) denote the free quark 
propagators with off-shell quark momentum p. In 
Fig. y(a) we plot the mixing coefficients 21, . . . , 24 
as a function of p^a^, where the vertical line cor- 
responds to about 2GeV. As mentioned in the 
previous subsection, this method can work under 
the condition Aqcd <S^ p <S^ a^^ ■ A strong p de- 
pendence of 2i, . . . , 24 below p « 2GeV indicates 
that non-perturbative contamination plagues the 
mixing coefficients up to around 2GeV. The Gri- 
bov uncertainty in Landau gauge fixing is another 
source of systematic errors. 

In the third method the mixing coefficients 
2i , . . . , 24 are determined to satisfy the contin- 
uum Ward identity up to 0{a) for the exter- 
nal off-shell quarks with momentum p pl| , p2[ : 
{V^,A^6LL[ssdd]ip)) = 2mq{PdLL[ssdd]{p)) + 
{i5A{OLL[ssdd]{p)}) +0{a). It should be stressed 
that this method can fix the mixing coeffi- 
cients without 0{mq) ambiguities for the (27, 1)- 
operator in the flavor 811(3)^x811(3)^, represen- 
tation which is relevant to Bk and K -^ Tnr[I — 
2] amplitude, while some ambiguities of 0{mq) 
remain in the case of the (8, l)-operator associ- 
ated with is: ^ 7r7r[/ = 0] decay ||. Figure |(b) 
shows the momentum dependence of the mixing 
coefficients. The Ward identity method requires 
the condition p ^ a^^ to avoid the cut-off ef- 
fects. We observe only weak dependence over 
a wide range p^a?<1.0, which contrasts to the 
NPMOM results in Fig. ||(a). In principle, how- 
ever, both results are expected to be consistent at 
the limit of large external quark momentum [p3| . 
The Gribov uncertainties in this method could be 
avoided by employing the Schrodinger functional 
technique. 

In Fig. |3| we summarize recent quenched results 
for Bk with the Wilson and SW quark actions. It 
is encouraging to be able to deal with the mixing 
problem without ChPT; still, large errors hinder 
the direct comparison with the KS result in the 
continuum limit. 



3. K ^ TTTT decay amplitudes 

There are two major difficulties in the lattice 
QCD calculation of X — > tttt decay amplitudes. 
The first one is the phase shift in tt-tt final state 
interactions, which cannot be directly calculated 
with real Green functions on a Euclidean lat- 
ticeQ. Secondly, we have to control the opera- 
tor mixing. For the Wilson quark action we have 

0cont ^ ^j^Cllatt _^ ^/^0/^ j^ (^^ _ mu){z3Sd + 

z'^inis -md)sj5d + Z5gosa^i,Ffj,i,d}, where O'^ are 
dimension-six operators arising from the explicit 
chiral symmetry breaking in the action. The mix- 
ing with sd, sj^d and g^isa^vF^^d is generated by 
the "eye" diagram, in which two quarks in 0± 
generate a quark loop by contracting with each 
other. The factors {nic — mu) and (m^ — md) 
are required by GIM mechanism and GPS sym- 
metry [g^ respectively. 

The simplest solution to overcome these dif- 
ficulties is calculating the amplitude {-n{p — 
Q)tt{p = 0)|O5?"t|K(pif = 0)) with rus = ma^. 
The final state pions at rest relative to each other 
generate no phase shift; the operator mixing is 
completely avoided by GPS symmetry and par- 
ity conservation. This method, however, requires 
GhPT to convert the amplitudes with the unphys- 
ical condition to those with the physical one. The 
next simplest method may be evaluation of the 
K ^ t: matrix elements with the KS quark ac- 
tion|33,M. Although the mixing with Oj_ can be 
avoided by the Ua(1) chiral symmetry retained in 
the action, subtraction of the lower dimensional 
operator is still necessary. This method also has 
to rely on GhPT. A new idea to avoid the operator 
mixing problem was recently proposed pq], which 
is an attempt to directly calculate the hadronic 
matrix element of the T-product of weak currents 
T[J^l{x)J'Il{^)] in the region a < |a;| < Aqcd- 
However, its practical feasibility is questionable 
even with 1 Teraflops machine because of the re- 
quirement of large cut-off. 

In recent years Golterman and his collabora- 
tors [Q,^ have studied the K -^ -kt: decay am- 
plitudes using GhPT and the quenched formula- 
tion of GhPT(QGhPT)|3| up to one-loop level, 
which enable us to investigate three systematic ef- 
fects contaminating the lattice QGD calculation: 
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Figure 5. K ^ tttt decay amplitude in A/ = 3/2 
channel with (filled symbols) and without (open 
symbols) one-loop corrections of QGhPT. Data 
are plotted as a function of lattice meson mass. 

unphysical kinematics, finite volume effects, and 
quenching. Only the first one can be estimated 
at the tree level. Although an 0(m|-/A^) un- 
certainty still remains in this calculation due to 
the unavailability of the 0{p'^) low-energy con- 
stant terms, we should note two points: First, 
finite volume effects can arise only from the loop 
diagrams. Second, we can investigate the differ- 
ence of the chiral properties between full QGD 
and quenched QGD, which emerge at the one- 
loop level. 

In Fig. we plot the physical X ^- 7r7r[/ = 2] 
amplitude as a function of the degenerate K and 
TT meson mass on the quenched lattice. They are 
obtained from the matrix element (7r(p = 0)7r(p = 
O)|0co„t|^(-^^ = 0)) calculated with the Wil- 
son quark action En,Bl|, where O™"^* is renormal- 
ized in NDR scheme at the scale 2GeV employing 
the tadpole improvement. The conversion to the 
physical matrix element is made by 



(^^|0+|X)phys 



2AB 



Y-{TnT\0+\K)utt,i^) 



where the factor Y denotes the one-loop contri- 
bution given by 



1 



Y 



{i^ur 



-104.73 -29.571og( 



^ + T£k7i-^^og{^)+F{AUL)] 



(4) 



F(Mt^L) represents finite size effects, whose 



expression is F{M^L) = 17.827/{M^L) + 
127r^/(M^L)^. The tree-level results(open sym- 
bols), which are obtained by setting m-^ — 
136MeV, ruK = 497MeV and F = 1 in eq.(|), 
are about two times larger than the experimen- 
tal value. However, once we include the one-loop 
effects with Y of eq.(Q) employing f^ = F^^ = 
132MeV and A = IGcV, the lattice results are 
reduced and are much closer to the experimen- 
tal value. The one-loop corrections are rather 
substantial. For example, the corrections are 
—40% for unphysical kinematics, —30% for finite 
volume effects, and -1-30% for quenching around 
Afjr w rriK at /3 — 6.1 with the 24^ spatial vol- 
ume(squares). It is also noted that the size depen- 
dence observed in the tree-level results is almost 
explained through the one-loop calculation. 

For the A/ = 1/2 amplitude, the one-loop cal- 
culation with QChPTB finds that the quenched 
matrix element {tt{p = 0)7r(p = 0)|O5?"*|/s:(pV = 
0)) has finite volume contributions in proportion 
to Mt^L, which arise from the rescattering dia- 
gram of the two pions. Although this fact could 
jeopardize the numerical calculation of (7r(p = 
0)7r(p = 0)\Of''*\K{pK = 0)) in the quenched 
approximation, fortunately, we have had no work 
so far. This dangerous finite size effect caused by 
the TT-TT interactions can be avoided by employing 
the A' — > TT matrix element. Its numerical cal- 
culation was tried by the OSU group employing 
the KS quark action[|2|. Although their result 
shows a consistency with the experimental value 
around a~^ « 2.8GeV, it is not convincing be- 
cause of their use of the tree-level ChPT and the 
statistical error of 40%. 

4. Proton decay amplitudes 

One of the most important features of the 
study of the baryon number violating processes 
is that the low energy effective theory is de- 
scribed in terms of SU(3)xSU(2)xU(l) symme- 
try, which enables us to make a model indepen- 
dent analysis. All dimension-six operators as- 
sociated with baryon number violating processes 
are classified in the four types under the require- 
ment of SU(3)xSU(2)xU(l) invariance p|,||. 
In the two-component notation of ref. ||44[| they 
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Figure 6. q^a^ dependence of the relevant form 
factor (circles) in {TT'^\{udii)uL\p), together with 
the form factor (triangles) employed in ref. M for 
comparison. Solid line denotes the fitting result. 

are given by o'^^^^^ = {d'^Ru'^R){q]aJjdL)ea/3je,j, 



O^aL = {€aL<l%L){< 



ldR)^al3j' 



Ri'dR 



V> 



(3) 
abed 

(4) 
ibcd 



i<lu^L(ljbL)i<lkcJndL)eaf3^e^nejk and O 

idaR'^bR)i''^lR^dR)£af3-f, whcrc a,(3,j are SU(3) 
color indices, i,j,k,n are SU(2) indices, a,b,c,d 
are generation indices and L and R de- 
note left- and right-handed fields. Specify- 
ing the decay processes of interest, in our 
case (proton, neutron)^ (Tr,/^) +{v, e+, /i^), we 
can list the complete set of independent ma- 
trix elements in QCD with the assumption 
of isospin symmetry: (7r°|eQ,^-y(u"(i^ ^)u]^|p), 
(7r+|e„^^(w"4jdlb), {K%^p^{u'^s%^)ul\p), 
{K+\e^p^{u^d%^)sl\p), {K+\e^p^{u"s%^)dl\p) 

and {K^\eai3'y{d°-s^^^)ul\p). All we have to cal- 
culate in lattice QCD are these matrix elements. 
Under the requirement of Lorentz invariance, 
we find that the matrix elements of O^ between 
the nucleon(iV) and the pseudoscalar(PS) meson 
can have two form factors: 

{PS{p)\0l\N\k)) = PL{W^{q^)+W,{q^)i)u^ , (5) 

where u^ denotes the Dirac spinor with spin state 
s. In the lattice calculation, fc = is chosen for 
the nucleon spatial momentum and p = k — q ^ 
for the PS meson. Although the Wq term in eq.(||) 
is irrelevant in the physical decay amplitude, be- 
cause its contribution is 0{mi) after the multi- 
plication of anti-lepton spinor, we have to disen- 
tangle these two form factors in the lattice QCD 
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Figure 7. Comparison of relevant form factors 
(circle) with predictions of tree-level ChPT (tri- 
angles). For the previous ChPT results (crosses) 
see text. 

calculation. As pointed out by JLQCDfl^], the 
previous papers IqM gave wrong results without 
realizing the existence of Wq term: they eventu- 
ally evaluated Wq + (pWq instead of Wq. 

In Fig. ^ the quenched results for WqW') in 
(7rO|(udfl)uL|p) calculated by JLQCdQ with 
the Wilson quark action at /3 = 6.0 are compared 
with Wo-|-9°Wq obtained by following the method 
in ref. g. The magnitude of Wo((3'^) is more than 
two times larger than that of Wq -I- (pWq. The 
value at q^ — O(open circle) is obtained by fit- 
ting the data with the function of cq -f CiQ^ -f 
c-i q'^ + c^niy^ci + c^rris . In Fig. we compare the re- 
sults obtained by the direct method(circles) with 
those by the indirect one(triangles) using tree- 
level ChPT, where the so-called a, (3 parameters 
are |a| = 0.015(l)GeV^ and \f3\ = 0.014(l)GeV^ 
We observe that both results are roughly com- 
parable, which could leads to the conclusion that 
the large discrepancy between two methods found 
in refs.p,p[ is mainly due to the neglect of the 
Wq{q'^) term in eq.(o). Another important point 
in Fig. Yj is that the JLQCD results are much 
larger than the tree-level ChPT predictions with 
|a| = \P\ = O.OOSGeV^ (crosses), which is the 
smallest value among various model estimations. 

5. Domain wall QCD 

If we can make the length of the fifth dimension 
N sufficiently large in the practical implementa- 
tion of the domain wall QCD on the lattice, the 



desirable features would emerge as: (i) no fine 
tuning for the chiral limit up to 0(l/ae~"^), 
(ii) 0{a'^) scaling violation up to 0{ae~°' ^) and 
(iii) no mixing for four-fermi operators up to 



0{g' 



-a"N 



). The flavor symmetry is retained at 



any number N. The KS quark action, however, 
can realize the above features without corrections 
except the flavor symmetry. The cost of numer- 
ical studies with the domain wall quark are the 
fifth dimension with the extension TV and choos- 
ing the optimal value for the domain wall height 
M. The RIKEN-BNL-Columbia group presented 
numerical studies for weak matrix elements rele- 
vant to K^-K^ mixing, K ^ tttt decays and the 



prediction of e' /e\2A\, assuming the naive chiral 
properties realized at the infinite iV|24|. In Fig. y 
we already plotted the quenched results for Bk 
by the RIKEN-BNL-Columbia group. Their re- 
sults are considerably smaller than the KS results 
at finite lattice spacings, and the discrepancy is 
likely to survive even in the continuum limit. At 
present it is not clear what causes this difference. 
Perturbative calculation in domain wall QCD, 
which was initiated in ref. [ [45| , has been applied to 
evaluation of the renormalization factors for the 
quark mass[^,^ and the bilinear operators |4q] 
consisting of four-dimensional quarks which were 
referred to as q{x) in ref. [|l2|. The Tsukuba 
group iQ shows that Zs = Zp = Z^^ and 
Zy — Za at N —t oo, which is expected when 
the chiral Ward-Takahashi identity hold exactly. 
The peculiar feature in the renormalization of 
the domain wall QCD is an appearance of the 
overlap factor (1 — |1 — M\^)Z^ for the four- 
dimensional quark fields. The one-loop coeffi- 
cient of Z^ is of 0(10 - 10^) for |1 - M|>0.1 
without mean-field improvement. This could be 
dangerous because current numerical studies em- 
ploy A/>1.6[p3. In the mean-field analysis, how- 
ever, they find that M should be replaced by 
M ^ M ^ 4(m — 1), where u"* corresponds to 
the expectation value of the plaquette, and the 
magnitude of the one-loop coefficient for Z^ is re- 
duced to 0(1) for |1-M|<0.8. This fact suggests 
that the mean field improvement is indispensable 
for the perturbative renormalization factors in do- 
main wall QCD. The p meson decay constant /p 
obtained with V^ = q(x)"ffj_q(x) and with the five- 
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Figure 8. fp^ and /^r obtained from the five- 
dimensional conserved current (circles) and the 
four-dimensional local current with (filled trian- 
gles) and without (open triangles) perturbative 
renormalization factors. 

dimensional conserved current V^°'^ gives a test- 
ing ground for the validity of perturbation the- 
ory ESI. Both results at f3 = 6.0 are compared in 
Fig.lsfa). They show a consistency within error 
bars once the perturbative corrections are applied 
to y. We also observe a similar situation for 
the case of the pion decay constant in Fig. @(b) . 
This is an encouraging result; still, the scaling vi- 
olation effects should be checked. The Tsukuba 
group also calculated the renormalization factors 
for the three- and four-quark operators associ- 
ated with proton decay and weak interactions ||49| . 
They find that the operators in domain wall QCD 
at infinite TV can be renormalized without any 
operator mixing between different chiralities, as 
opposed to the Wilson case. 

The good chiral properties of the four- 
dimensional quarks in domain wall QCD observed 
in the numerical simulations and the perturbative 
calculations can be understood from the point of 
view of the exact chiral symmetry based on the 
Ginsparg- Wilson relation. The partition function 
of domain wall fermion with the subtraction of 
the Pauli-Villars field at finite N can be written 
as a determinant of the truncated overlap Dirac 
operator _DArpffl which satisfies the Ginsparg- 
Wilson relation in the limit of infinite N. The 
authors of ref. p^ find that the effective Dirac 
operator D'^ for the four-dimensional quarks in 
domain wall QCD with finite N, which is ob- 
tained by integrating out the A^— 1 heavy fermion 
fields, is related to Dn: Z3^ {x,y) + S{x,y) ~ 



DN^^{x,y). With the use of this relation they 
show that the Green functions consisting of q{x) 
and q{x) in domain wall QCD with finite N can 
be related to those consisting of the fermion fields 
described by the truncated overlap Dirac opera- 
tor. 

6. Conclusions 

The sizable magnitude of the 0{a^) uncertain- 
ties found in the systematic study of Bk with 
the KS quark action makes manifest the neces- 
sity of using the non-perturbative renormaliza- 
tion method. The essential problem in the non- 
perturbative renormalization of QCD, which is 
the question how to relate the high-energy pertur- 
bative regime to the low-energy hadronic scale, is 
overcome by the Schrodinger functional method. 
For the Wilson Bk the operator mixing problem 
is now managed non-perturbatively without the 
use of any effective theories; while we still need 
to reduce statistical errors. 

In K -^ TTvr studies one-loop corrections in 
QChPT can almost explain the discrepancy be- 
tween quenched lattice results and the experimen- 
tal value for the A/ — 3/2 transition. 

Essential progress has been made in evaluating 
the proton decay amplitudes. The first correct 
calculation shows that the values of the ampli- 
tudes are comparable with the tree-level ChPT 
predictions. 

The application of the domain wall quark for- 
mulation to the calculation of the weak matrix 
elements is a fascinating issue. Pilot numeri- 
cal studies have been pursued and perturbative 
renormalization factors for the relevant operators 
are now available. 
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